Asymptotically noise decoupling for Markovian open quantum systems 
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The noise decoupling problem is investigated for general N-level Markovian open quantum sys- 
tems. Firstly, the concept of Cartan decomposition of the Lie algebra su(N) is introduced as a 
tool of designing control Hamiltonians. Next, under certain assumptions, it is shown that a part of 
variables of the coherence vector of the system density matrix can be asymptotically decoupled from 
the environmental noises. The resulting noise decoupling scheme is applied to one-qubit, qutrit and 
two-qubit quantum systems, by which the coherence evolution of the one-qubit and qutrit systems 
can always be asymptotically preserved, while, for two-qubit systems, our findings indicate that 
evolution of some variables can be preserved only for some initial states. 
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I. INTRODUCTION 

In recent years, quantum information science 
hasElIlllHlIlBll been a growing field which inter- 
ests many researchers for potential high speed quantum 
computation and high security quantum communication. 
Decoherence is commonly recognized as the main bot- 
tleneck. Various schemes have been proposed to reduce 
such unexpected effects. In principle, there are two 
classes of schemes — open loop and closed loop strategies 
depending on the use of measurement and feedback. 
Open loop strategies include quantum error-avoiding 
codesBlia, Bang-Bang control^ El El El, open 
loop optimal controlEl El El anc j open loop coherent 
controlE^ El ; while closed loop strategies include quan- 
tum error-correction codes^'^E'^E'^l anc j quantum 
feedback control^H!. 

Though various strategies have been proposed, none 
of these strategies is satisfying to suppress decoherence 
for iV-level Markovian open quantum systems: quantum 
error-correction codes and error-avoiding codes use sev- 
eral physical qubits to encode one logical qubit, which is 
too luxurious under existing conditions; BangBang con- 
trol strategy is inapplicable in the fully Markovian regime 
as pointed out by LidarEl ; open-loop optimal control 
strategy can only partially decouple quantum systems 
from the environmental noises; quantum feedback con- 
trol strategy requires complex feedback control apparatus 
and it is valid only for special physical systems. Thus, for 
the decoherence suppression problem of iV-level Marko- 
vian open quantum systems, more system analysis and 
control methods should be introduced. 

The closest work to ours can be found in Ref. 19], 
where open loop coherent control is applied to decoher- 
ence suppression for single-qubit Markovian systems. It 



is shown that the x-axis and y-axis variables of the Bloch 
vector can be exactly decoupled from the environmental 
noises, which means that the coherence of the quantum 
state can be well preserved. However, this scheme re- 
quires solving a time-variant linear ordinary differential 
equation (ODE) to obtain the open loop control laws, 
by which analytic control laws can be obtained only for 
phase damping decoherence under a strong assumption 
that the x-axis and y-axis variables of the Bloch vector 
keep constant. For more general cases, only numerical 
control laws can be obtained. Furthermore, divergence 
of the control fields may occur in this strategy as pointed 
out by the authorsEl . In this paper, we propose a more 
general noise decoupling strategy for TV-level Markovian 
open quantum systems based on the Cartan decomposi- 
tion of the Lie algebra su(N). Open loop controls are 
designed to asymptotically decouple the state variables 
from the environmental noises. The strategy loses some 
precision but is easier to be fulfilled. 

The paper is organized as follows. In section|TTl Marko- 
vian open quantum systems are formulated in the coher- 
ence vector representation and the concept of the Cartan 
decomposition of su(N), together with three important 
assumptions, is introduced. Our main results are pre- 
sented in section IIII1 where the scheme are also applied 
to one-qubit, qutrit and two-qubit systems. Further dis- 
cussion and conclusion are drawn in section ITVl 
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II. PRELIMINARIES 

Consider an N-level Markovian open quantum control 
system in the following master equation form: 

n m 
i=l j=l 

p{to)=Po, (1) 
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where the Planck constant h is assigned to be 1; p refers 
to the system density matrix; H Q and iij, i = 1, • • • , n 
represent, respectively, the free Hamiltonian and the con- 
trol Hamiltonians adjusted by the control parameters 
m, i = 1, • • • , n. The Lindblad super operators 

B[Lj]p = Lj P L\ - \l]l ]P - \pL\L h 

characterize the damping channels and the positive con- 
stants Tj denote the damping rates of the corresponding 
channels. 

The differential equation ([1]) is actually a complex ma- 
trix differential equation which is hard to be analyzed. 
Therefore, we will convert it into a real vector differen- 
tial equation. For this purpose, an orthonormal basis 
{Oo = ~^fil> fy}j=i,— ,N 2 -l with respect to the matrix 
inner product (X,Y) = tr(X^Y) should be introduced 
first, where I is the N x N identity matrix and Q'jS are 
N x N Hermitian traceless matrices. The system density 
matrix p can then be expressed as: 

1 W2_1 1 
p = — I + m ^ ~jjl + m-n, (2) 

i=l 

where m e R*"- 1 is the so-called coherence 
vector^ EE El [IE [III, OH] ^ ^ ^ ^ g cage ^ ^ q uan ^ um 

control system (fTJ) can be reexpressed by a differential 

equation on M. N _1 : 

n 

fn(t) = O m(t) + ^ UiOim(t) + Dm(t) + g, 

i=l 

m(to) = mo, (3) 

where Oo,Oi € so(N 2 — 1) are, respectively, the adjoint 
representation matrices of —iHo, — ii^OH and too is the 
coherence vector of po- The term 11 Dm + g" comes from 
the decohering process represented by the Lindblad terms 
in |T]) in which the dissipative matrix D is semi-negative 
denned, i.e., D < 0. 

The above approach is a generalization of the well- 
known Bloch vector representation for two-level quantum 
systems. Physically, the length of the coherence vector 
represents the amount of coherence in the quantum state. 
Concretely, it is unit for the pure state and shorter for 
the mixed state. Notice that the control ^ UiHi drives 
the quantum state along a sphere on which coherence is 
conserved, while the decohering operators pull the vector 
towards the equilibrium state, e.g. the ground state in a 
spontaneous emission. 

To simplify the equation © and facilitate our discus- 
sions, it is useful to discuss the choice of the matrix basis 
{-7=1, ilj}j =1 ... N 2_ 1 . In this regard, we introduce the 

so-called Cartan decomposition^ 03 03 of the Lie al- 
gebra su(N) as follows: 

su{N) = p © e, [e, e] C e, \p,p] C e, \p,e]c p. (4) 



Notice that su(N) is a N 2 — 1 dimensional Lie al- 
gebra of all traceless skew-Hermitian N x N matrices, 
hence the basis matrices of p and e can be expressed 
as {-if2f}i=i,..., m and {-iflf}/= m +i ) ... ,jv 2 -i where fi? 
and £lf are traceless Hermitian matrices. Further, ev- 
ery traceless Hermitian matrix is a linear combination of 
, f^}. Correspondingly, the density matrix p can be 
represented as: 

1 rn N 2 -l 

p=N I +Y, m i n i+ E ( 5 ) 

i=l l=m+l 

As will be shown later, the matrices {^f}j=i,... ,m play 
central roles in our control strategy. In fact, we will 
choose the control Hamiltonians Hi from these 
matrices. 

It should be pointed out that the Cartan decomposi- 
tion always exists. In fact, we can obtain a trivial de- 
composition if we let e = su(N). The decomposition is 
also not unique. 

Before proceeding our discussions, we introduce three 
assumptions for the equation Q: 

(HI) Complete decoherence condition: 
[O ,D]=0, O g = 0: 

(H2) Convergence condition: D < 0; 
(H3) -iH e e. 

(HI) has an important physical interpretation's that 
the stationary distribution p^ of the uncontrolled system 
satisfies [HojPoo] = (see the appendix of Ref. [17| for 
a rigorous proof). In other words, in the energy repre- 
sentation, the off-diagonal entries of the stationary sys- 
tem density matrix p^ disappear as a result of deco- 
herence. (H2) is introduced to guarantee the existence 
of the convergent solution of (J3J). In fact, we can al- 
ways make D < with the aid of the feedback control 
modificationE3 . (H3) is easy to be satisfied under special 
choices of the Cartan decomposition. 

III. ASYMPTOTICALLY NOISE DECOUPLING 
STRATEGY 

The term 11 Dm +g" in equation ([3]) is the environment- 
induced dissipative term which destroys coherence of the 
quantum states. Our target is to select the control Hamil- 
tonians Hi and design the corresponding controls m to 
force the trajectory m(t) of the equation ([3]) as close as 
possible to the target trajectory m°(t) which is the solu- 
tion of the following unperturbed system: 

m(t) = Oom(t), m(to) = mo- (6) 

It has been demonstrated that the Markovian open 
control system ([3]) is always uncontrollable 1^3 . There- 
fore, one can never track the target trajectory precisely. 
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However, we will show that certain variables of m(t) 
may asymptotically tend to the corresponding variables 
of m°(t) by properly designed control laws. That is to 
say, these variables can be asymptotically decoupled from 
the environmental noises. 

In fact, according to ((5]), the coherence vector m can 
be divided into two parts: to = (to 1 , m 2 ) T , where 



! '"-ml i 



— (m^ l+1 ,- ■ ■ ,m £ N 2_ 1 ) 1 



From lemma lA~2| we have Oq = diag^l 1 ,0^ 2 ), where 
Oq 1 , Oq 2 are respectively, to and N 2 — m — 1 dimensional 
square anti-symmetric matrices, then the target trajec- 
tory can be written as: 



where 



m°(t) = e 0o(t -* o) TO = (m w (t),m 20 (t)) T , 

m w (t) = e^-^ml, 
m 20 (t) = e^^-^m 2 , 

and too = (toJ,TOq) t . In this case, the vector m 1 ^) can 
be driven to the corresponding target trajectory m w (t). 
In fact, we have the following theorem: 

Theorem 1 Suppose the assumptions (HI), (H2) and 
(H3) are satisfied and the control Hamiltonians Hi in 
(QJ) are the basis matrices {Of} corresponding to p. The 
following control law 



(«i,-,«m) r = e-°o 11 ( t - t0 ^ 



(7) 



steers the control trajectory m(t) of the equation fj) 
asymptotically to the stationary solution: 

m°°(t) = (e o "(*- tD >mS,e o "< t - t0 >»7) T , (8) 

where the constant vectors £ £ R m and rj £ M. N ~ m ~ 1 are 
the solutions of the following nonlinear algebraic equa- 
tion: 

m 



i=l 



»=i 



where 



D 



+D22V + .92 = 0, 



D21 D22 



(9) 
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9i 

.92 



Rigorous proof of the theorem is left to appendix 
Theorem[T]shows that the state variables of the coherence 
vector of p corresponding to p can be asymptotically de- 
coupled from the environmental noises. Note that Oq 1 is 
an TO-dimensional anti-symmetric matrix which has only 



zero or pure imaginary eigenvalues, the control laws 
presented in ([7]) are in fact sinusoidal signals which 
are usual for electro-magnetic fields used in the labora- 
tory. 

Next, we will apply our control strategy to three typ- 
ical systems in quantum information science: one-qubit, 
qutrit and two-qubit systems. 



A. One-qubit quantum systems 

Consider the one-qubit systems which are fundamen- 
tal in the quantum information science. Assume that the 
free Hamiltonian Hq = ljo~ z where o~ z is the z-axis Pauli 
matrix, then the Cartan decomposition of the corre- 
sponding Lie algebra su(2) can be chosen as su{2) = p(Be, 
where the basis of p and e are {— i-^=a x , — i-^=a y } and 

{— i-j=a z } respectively. 

Now, for the coherence vector representation of p, i.e., 



y V2 



m z7=2 * z 



we can conclude that the variables m. x = -j=tr(o~ x p) and 
m y — -j=tr(a y p) can be asymptotically decoupled from 
the environmental noises if the assumptions (HI), (H2) 
and (H3) are satisfied. 

As an example, we study the one-qubit amplitude 
damping decoherence modelSI which can be used to de- 
scribe spontaneous emissions of the two-level atoms. In 
this case, the control system can be expressed as the fol- 
lowing master equation: 

p = -i[uJo- z + u x a x + UyO-y, p] + TD[o--]p, (10) 



where the controls and lu £ 



denotes the 
o~ x — ia y is the lowering operator of 



Rabi frequency; ti- 
the two-level system and T > denotes the decoherence 
rate, e.g., the damping rate of the spontaneous emission 
process. To simplify the calculations, we use the well- 
known Bloch vector representation: 

p = -(I + m x a x + m y <j y + m z a z ) : 

where m x = tr{a x p), m y — tr{a y p) and m z = tr(a z p). 
Note that this representation is different from the co- 
herence vector representation only by a trivial constant 
multiplicative factor. 

The master equation pH|) can be converted into the 
following equation in the Bloch vector representation: 

to = u)O z m + u x O x m + UyOym + Dm + g, (11) 

where to = (m x , m y , m z ) T £ K 3 and 







1 



-1 







4 



-1 



O z 







are the basis elements of the Lie algebra so(3) of the 3 
dimensional orthogonal group; D = diag(—^k,—i>,—T) 
and g = (0, 0, — T) T come from the decoherence process. 
It can be verified from (flQ|) and (fTTjl that the assumptions 
(HI), (H2) and (H3) are all satisfied. 

From the equation [O ,Of] = J2j(Oo 1 )ijOj and 

[0 Z ,0 X ]=Oy, [0 Z ,Oy}=-0 X , 



it can be easily computed that Oq 1 



Ac- 



cording to theorem [1] the controls can be designed as: 



e -o^(t-t )^ 



cosuj(t — to)£i — sinw(t 
smuj(t — to)£i + cosw(t 

Acos(uj(t — to) + <f. 
A sin(aj(t — to) + <fi 



*o)6 



under which m x . 
of the unperturbed system m 
£ = (£i, £,2) T can be solved by 



tend to the corresponding variables 



uiO z m. The constant 
) as follows: 



Tm 0y 

~w 

Tm 0l 
2C 2 



1±a/1 



2C 2 



1 ± 



1 - 2C 2 



where mo 

0) = m 0x 



(mo x ,m y,mo z ) 



is the initial state and 



l 0y 



represents the initial coherence in the 
quantum system. It can be observed that Cg should be 
no larger than ^ to guarantee the existence of the solution 
of |9]) . Amplitude and phase of the control fields can then 
be calculated as: 



e 2 = — 



r 

2Co 



arctg— = arctg 



l ±v /l- 

m 0x \ 
m y) 



2CI 



Let u> = 3/tq, to — and the initial state 

1 T V2 V2 
A) = — + 

where To is a time constant which is introduced to obtain 
dimensionless evolution time. For concrete systems, To 
can be determined by the system time scale, e.g. relax- 
ing time of systems. With simple calculations, it can be 



shown that u x = sin(3t/r ), u y = 
Simulation results of the variables m x 



U x , U 



,V5 



rcos(3t/T ). 
and controls 



are shown in Figure [T] and [21 
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FIG. 1: Temporal evolution of (a) m x and (b) m y : the asterisk line 
represents the controlled trajectory; the plus-sign line is the uncon- 
trolled trajectory; the solid line is the target trajectory. 




FIG. 2: Temporal evolution of controls: the solid line is u x and the 
plus-sign line is u y . Here, to obtain dimensionless quantities, u x and 
u y are divided by the decoherence intensity T . 



Since the coherence C 2 = m 2 , + m 2 in p is determined 
by m x and m y , the coherence of the state will vanish 
completely without control. The figure shows that the 
controlled trajectory tends asymptotically to the target 
trajectory where the two trajectories are so close that 
they almost coincide together, which implies that the co- 
herence of the state is asymptotically preserved with our 
control strategy. 




FIG. 3: Temporal evolution of controls: the solid line is u x and the 
plus-sign line is u y . Here, to obtain dimensionless quantities, u x and 
u y are divided by the decoherence intensity T . 



In the study of Lidar et al on the same problem in Ref. 
[HI, the state variables m x and m y are exactly decou- 
pled from the environmental noises under the following 
feedback-like control laws: 



-m x , 



(12) 



■1,„ : tf ' 2m. 
which are to be substituted into Ijlip to get explicit open 



5 



loop control laws. This leads to the following equations: 



itt-x (t) — (t) — mox cos Lot + may sin tot, 
m y (t) = rriy(t) = —mo x sinujt + mo y cos art, 
m z rn z = —Tm 2 — Tm z 



r r< 2 



(13) 



l 0y 



where m = (m 0x ,moy,m 0z ) T and Cq = w,q x 
Generally speaking, the last equation of (fT3"| has no an- 
alytic solutions and we can only obtain numerical solu- 
tions. With the same parameters in the above example, 
we can obtain plots of controls u Xl u y in Figure [3J 

Compare Figure [2] and [3j it can be shown that the 
control laws in Figure [3] are more complex, and they ap- 
proach to our laws asymptotically. Furthermore, a di- 
vergent solution of the last equation of (fT3"|) may lead 
to divergent control fields. With simple calculations, it 
can be shown that solution of the last equation of (fT"3"]) is 
convergent if and only if 



1 



< 



Figure 2] shows the divergent control fields when the ini- 
tial state is chosen as m — (^, 0, -^) T . 



the environmental noises. The plot of the difference be- 
tween the trajectories obtained by our strategy and Li- 
dar's strategy (ideal trajectories) are shown in Figure [5] 



B. Qutrit quantum systems 

Consider the three-level atoms^GElGSISl with V con- 
figuration as shown in Figure [6j where |0), |1) and \2) 




FIG. 6: Three-level atoms with V configuration 
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Divergent control fields for the initial state 



FIG. 

to obtain dimcnsionlcss quantities, 
hcrcncc intensity T. 



The solid line is u x 



and the plus-sign line is u y . Here, 
and u y are divided by the deco- 



are the eigenstates of the free Hamiltonian of three-level 
atoms with eigenvalues Eq < E± < E2 respectively. The 
two excited states |1) and |2) are coupled to the ground 
state |0). 

The total Hamiltonian H is expressed as H = Ho + Hd 
where Ho is the free Hamiltonian and Hd is the control 
Hamiltonian. The free Hamiltonian Ho can be written 
as: 

^o = E |0)(0| + £ 1 |l)(l|+ J B 2 |2)(2|. 
The control Hamiltonian 



A m 




H d = ffio|l><0| +^ |0)(1| +92o|2)(0| + <? 2 *o|0><2| 

represents the interaction between the atoms and the 
driving electromagnetic fields. The complex coefficients 
310)320 can be adjusted by the amplitudes and phases of 
the driving fields. 

Consider the open system in which only amplitude 
damping decoherence channels induced by spontaneous 
emissions are introduced. The master equation model is 
expressed as: 



FIG. 5: Plot of Am,, = y(mj — m^) 2 + (mj — m y ) 2 where m], 
and m y are respectively x, y entries of the trajectory obtained by our 
strategy and m 2 , and m 2 arc x, y entries of the trajectory obtained by 
Lidar's. 

It should be pointed out that, though the control laws 
obtained by our strategy are simple and the divergence 
problem does not occur, our control laws can only asymp- 
totically, not exactly, decouple the state variables from 



p = -i[H + H d , p) + r^KUp + T 2 D[ao 2 ]p, (14) 

where ov^ = \0)(i\ is the lowering operator from the ex- 
cited state \i) to the ground state |0). The two Lindblad 
terms TiD[<Joi\p represent the transition from \i) to |0) 
caused by the spontaneous emission process. 

For three-level systems, the matrix basis of the cor- 
responding Lie algebra su(3) = {—i£lk}k=i,— ,& can be 
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chosen as: 



. , o \ 1 ( 

o, =-= o o i , o 2 = 7 = -i 

V 2 \ o 1 / v 2 \ i 

. , \ , f 1 

n 3 = -= o i o ,n 4 = T i o o 

V 2 1 -1 v 2 I o 



1 



-i 



Q 5 = —= \ i 
v2 \ o 



1 

n 6 = I o o o 

V 2 \ 1 



1 



0-i 



f2 7 = —j= 
v 2 \ i 



1 



-2 



v 6 \ 1 



It can be verified that sn(3) = p (B e is a Cartan decom- 
position where 

p = {—{0,4, —iCls, — iQ,6, —1CI7}, 
e = {-iQi, -iQ 2) ~*^3, -i^s}- 

With the basis matrices {fife}, the free Hamiltonian 
can be rewritten as: 

it o i n Eo + Ei + E 2 
H Q = uj 3 il 3 + us^k H jj i, 

where 

^3 = ^{E X - E 2 ), u 8 = ^-(E 1 + E 2 - 2Eo). 
2 o 

Since the constant energy (Eq + E x + I?2)/3 only con- 
tributes a global phase to the system state, it is sufficient 
to consider the following traceless free Hamiltonian 

Furthermore, let g 10 = ^(u 4 + m 5 ), .920 = ^("6 + ^7)- 
The control Hamiltonian can be written as: 

H d = U4SI4 + + u 6 fl& + u 7 Q 7 , 

where u[s are the controls to be designed. 

Now, from the coherence vector representation of p, 
i.e., 

1 8 

p = -I + y^TOjOj, mi = tr(£lip), (15) 



we obtain the following coherence vector representation 
of the master equation (jT4j) : 



rh = lu(<jj?,Oz + aJsOs)m + UiOim + Dm + g, (16) 

i=4 

where = ad(-ifti). It can be verified from (TH)) 
and (frnj) that the assumptions (HI), (H2) and (H3) 
are all satisfied. According to theorem [I] the variables 



rrii, i = 4, 5, 6, 7 in the equation (|15p can be asymptoti- 
cally decoupled from the environmental noises under the 
following control law: 



(u 4 ,U 5 ,U 6 ,U 7 ) 



T _ p -O?(t-t ) 



where 



/ 



o, 



11 



w 3 



2 u 

#00 





V 



\ 




f / 



+W 8 



/ ^ 

2 u 






\ 







V 






& 

u 2 

vs. 

2 



/ 



and £ can be numerically solved from the equation ([5]). 
Unlike one-qubit quantum systems, the algebraic equa- 
tion © has no analytic solution in this case. 

Let r x = r 2 = r, t = o,e = -*m, El = __ug i 

E% 



tt^ and the initial state be the mixed state: 



PO 




where tq is a time constant which is introduced to obtain 
dimensionless evolution time. For concrete systems, r 
can be determined by the system time scale, e.g. re- 
laxing time of systems. With simple calculations, it 
can be shown that u 4 = — 0.7063r sin(4.4t/ro), 1*5 = 
-0.7063rcos(4.4i/r ), u e = -0.7063r sin(2.5111t/r ), 
W7 = — 0.7063rcos(2.5111</r ). Here, the amplitudes 
and initial phases of the control fields are obtained by 
numerically solving the equations (J5J). Simulation results 
of the variables rrij and controls Uj, i — 4, 5, 6, 7 are shown 
in Figure E| and [H] 



Recall that the coherence C§\ = tci\ 



between 



|0) and |1) is determined by 7714 and 7715; the coherence 
Cq 2 = mg + niy between |0) and |2) is determined by m 6 , 
7777. The coherence in p will vanish completely without 
control according to Figure [7] The controlled trajectory 
is driven so closely to the target trajectory that they al- 
most coincide together, i.e., the coherence between the 
excited state |z) and the ground state |0) can be asymp- 
totically preserved with our control strategy. 

Just like one-qubit systems, if we want to exactly 
decouple the corresponding variables from environmen- 
tal noises, only numerical control laws can be obtained. 
However, with the same parameters in the above exam- 
ple, we find that the corresponding time- variant ordinary 
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FIG. 8: (a) Plot of U4 and u$ where the solid line is U4 and the plus- 
sign line is u$; (b) Plot of uq and uy where the solid line is uq and the 
plus-sign line is 117. Here, to obtain dimcnsionlcss quantities, u\ are 
divided by the decoherence intensity F. 



differential equation to obtain the numerical control laws 
has no solutions. In fact, in order to solve the equation, 
we must calculate an inverse matrix and this matrix is 
singular for the initial state in the example. It means 
that exactly decoupling strategy fails for this example. 

Finally, the results obtained by our control strategy 
can be directly extended to the A-type three level atoms 
and other kinds of decoherence channels including phase 
damping decoherence channels and depolarizing decoher- 
ence channels, as long as the assumptions (HI), (H2) and 
(H3) are satisfied. 



C. Two-qubit quantum systems 

For two-qubit systems, the corresponding Lie algebra 
su(4) has a Cartan decomposition su(4) = p © e, where 

r 1 

P = Xi^i ® a j\hJ = x,y,z\, 

e = 7j<Ji®I\i,j = x,y,z}. 

Now the coherence vector representation of p can be 



written as: 

i 3 

where 

m\ = \tr((Ji <g> I)p, m 2 = \tr{I <g> aj)p, 
m\] = ^tr(ai®aj)p. 

It is shown that the two-qubit variables {mj 2 } can be 
asymptotically decoupled from the environmental noises 
if the assumptions (HI), (H2) and (H3) are all satisfied. 

As an example, consider the two-qubit independent 
amplitude damping decoherence model which describes 
two atoms that simultaneously undergo spontaneous 
emissions. The control system is described by the fol- 
lowing master equation: 

p = -i[H + ^2 UijHij , p] 

+T l D[^a 1 _®I]p + T 2 D[h®a 2 _] Pl (18) 

where H Q = u\\o\ ® I + <8> a 2 and H t j = \a\ ® 
a 2 , i,j = x,y,z, are the free and control Hamiltonians 
respectively; a % _ = a l x — ia y , i = 1, 2 is the lowering op- 
erator of the i th subsystem and the positive coefficients 
Tj denote the corresponding decoherence rates. The two 
Lindblad terms represent the amplitude damping deco- 
herence channels of the two subsystems. 

Rewriting the master equation (I18p in the coherence 
vector representation, one can verify that the system sat- 
isfies the assumptions (HI), (H2) and (H3). Therefore, 
our control strategy ((7|) can be applied. Unlike the one- 
qubit case, the quadratic algebraic equation @ does not 
have analytic solutions, however, we can obtain numeri- 
cal solutions of £. 

Let u>i = uj 2 = 1/t , t = 0, Ti = T 2 = T, where r 
is a time constant which is introduced to obtain dimen- 
sionless evolution time, and the initial state be the mixed 
state: 

Pa = 2 ■ \ I+ g l^oX^ol, 

where \<f>o) = -^(|00) + |11)) is the maximally entangled 
Bell state. In this case, among the two-qubit variables 
{mjf} in the equation (|T7|1. only m x 2 , m^ 2 ,, m x 2 x , m^ 2 
and mil are non-zero. The simulation results of these 
variables are shown in Figure [5J It is shown that the 
uncontrolled trajectories of roH, mz 2 ., m x 2 x and m^ 2 and 
m 12 z evolve away from anticipated values. With controls 
plugged in, the controlled trajectory tracks asymptoti- 
cally the target trajectory so close that they almost co- 
incide together. 
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FIG. 9: Temporal evolution of (a) m™, (b) m y \, (c) m**, (d) m y 2 y 
and (e) m^: the asterisk line denotes the controlled trajectory; the 
plus-sign line is the uncontrolled trajectory; the solid line is the target 
trajectory. 



Further, with simple calculations, it can be shown that 
onlyu^ = -ircos(2i/r ), u yx = -±r cos(2i/r ), u xx = 
— s'm(2t / To) , u yy — |T sin(2t/ro), and u zz — jT are 
non-zero. Plots of the controls are shown in Figure HOI 



See Figure [Til for an example. 
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FIG. 11: Temporal evolution of (a) m^, (b) mf x , (c) m^, (d) mf y 
and (e) m zz with wi - QJ2 - l/^Oj Ti — T2 — T, fo — and the 
initial state being the maximal entangled Bell state po — |0o)(0oli 
where |0o) — -^^= ^] 00} + The asterisk line denotes the controlled 

trajectory; the plus-sign line is the uncontrolled trajectory; the solid 
line is the target trajectory. 




IV. DISCUSSION 



FIG. 10: (a) Plot of u xy and u xx where the solid line is u xy and the 
plus-sign line is u xx ; (b) Plot of u yx and u yv where the solid line is u yx 
and the plus-sign line is u yy . Here, to obtain dimcnsionlcss quantities, 
Uij are divided by the dccohcrcnce intensity T. 

Like one-qubit and qutrit systems, if we want to 
exactly decouple the variables from the environmental 
noises, only numerical control laws can be obtained. 
With the same parameters in the above example, we find 
that the corresponding time- variant ordinary differential 
equation to obtain the numerical control laws has no solu- 
tions, which means that exactly decoupling strategy fails 
for this example. 

It should be pointed out that the algebraic equation (J9]) 
does not always have solutions. The existence depends 
on the initial state. In any case, we can always find the 
least-squared solution of the algebraic equation © by 
solving the following minimization problem: 

rmnF 1 (Z,r l ) 2 + F 2 (Z,r ) ) 2 , 

where Fi(£,rf), i — 1,2 are given in ([5]). With such con- 
trols, we can only partially recovers the target trajectory. 



In subsection IIII CI it has been pointed out that 
the two-qubit variables {mj 2 } can be asymptotically 



decoupled from the noises. In Ref. [37|, we pro- 
posed a multipartite mixed-state entanglement measure 
modified from Jaeger's Minkowskian norm entanglement 
measure^ Hl> S3, which, for two-qubit states, is defined 



E(p) = 



i 12 Y - - 



The entanglement measure is only related to the two- 
qubit variables {m\j}. For this reason, it is reasonable 
to expect that the entanglement of states can be asymp- 
totically preserved by our control strategy. 

Unfortunately, we find that, for most entangled states, 
the algebraic equation (J9j) has no solution and our control 
strategy can not preserve entanglement completely under 
the least-squared solution. It should be further studied 
to what extend our control strategy may help to preserve 
entanglement. 

Our asymptotical noise decoupling strategy applies 
control Hamiltonians from the Cartan decomposition of 
the Lie algebra su(N) to decouple the systems from 
the environmental noises under reasonable assumptions. 
Such control may not be applicable in laboratory in 
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present condition, especially for the two-qubit example, 
but it still provides useful hints in systematic design of 
decoherence control. 

The construction of the Cartan decomposition is 
essential in our scheme. Since the decomposition is not 
unique, the finding of a "good" decomposition to achieve 
the expected control performance is an interesting 
problem that needs further research. 
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APPENDIX A: PROOF OF THEOREM [JJ 

Before presenting the proof of the theorem, we first 
introduce two lemmas: 

Lemma A.l I^i/ Consider the following time-variant 
linear system: 

x{t) = A(t)x(t), x(t ) = x a . (Al) 

If there exist a N x N matrix Q{t) and positive numbers 
77, y,, v such that: 

rjl < Q(t) < nl, 
A T (t)Q{t) + Q(t)A(t) + Q{t) < -vl, 

we have: 

\x{t)\ 2 <^e-^-^\x \\ 
V 

The following lemma shows that, if the control Hamil- 
tonians Hi are chosen as fif in the control system 
([3]) has a simple structure. 

Lemma A. 2 Let Oq = ad(—iHo), where ad(A) is the 
adjoint representation matrix of A, then 

(1) we have Oq — diag{C<Q l , Oq 2 ), where Oj 1 , O 22 are 
respectively, m and N 2 — m—1 dimensional square 
anti- symmetric matrices. For 0\ — ad(-ifl^), we 
have: 

p _ ( O} 2 \ 

(2) we have the following equation: 

m 

e -O a {t-t Q ) qP e Oa(t-t ) _ y^^-Op 1 (t-t )j ,.QP ^ 
3=1 

where (e~°° ^* _t °^)ij is the ij th entry of the tri- 
dimensional matrix e~°« ( i— *°). 



Proof. 

(1) Corresponding to the Cartan decomposition Q, we 
have the following decomposition of the represen- 
tation space M. N ~ 1 : 

R^ 2 " 1 =]RP©]fr, (A2) 

where 

W ={(oi,.-- ,a m ,0,--- ,0)|oi eM}, 

M e = {(0,--- ,0,a m+1 ,--- .ajv.-iJlai GK}. 

From the assumption (H3), Oq can be written as: 

N 2 -l 

Oq = uiadi-mt). (A3) 

l=m+l 

According to the Cartan decomposition (TJ|, the 
special structure of Oq and Of can be easily veri- 
fied from the fact that ad(p) maps W into R e , W 
into W, and ad(e) maps W into R e , W into W. 

(2) From (|A"3j) and the equality ad([A,B\) = 
[ad(A), ad(B)}, it can be shown that: 

m 

[Oo,Of]=E(°o% , P - (A4) 

3=1 

In fact, it can be deduced that: 

N 2 -l 

[Oo,0?)= w,[ad(-<nf),ad(-*n?)] 

l=m+l 

N 2 -l N 2 -l m 

= Y, uiad([-iQi, — ifif]) = ^ y]wicnjad(—i£l%) 

l—m+l l—m+1 j — 1 

rn I N 2 -l \ m 

= E E ^ oj = E(°")«°J. 

3=1 \i=m+l / 3=1 

where {cijfc} is the structure coefficients of the Lie 
algebra su{N). From (|A4|1 . it can be easily verified 
by induction that 

m 

[(-Oo) (fe) ,Of]-E((- o 1 ) fc ) y ^, 

3=1 

where [A®,B] = [A, [A^ 1 ), 5]] and [A^\B] = B. 
Now, from the equality 

oo 

e A Be- A = Y±[A^,B], 

i=0 



10 



we have: 



e -Oa{t-t )QP e o (t-to) 
E^[(-Oo(t-to)) (fe) ,Of] 

k=0 
oo rn . 

EE^((- o n (*-*o))%Of 



fc=0 3=1 



= E( e ~°° 11( '~ to) )^- 



where mg° = (m,Q,ri) T . From the complete decoher- 
ence condition (HI), we have [e° ^ to \D] = and 
e o (t-t )g _ ^ w hich results in 



E "iO? e^-^m^ + e° o{t - to) {Dm™ + g) = 0. (A7) 



Substituting ([7]) into the first term in (|A7p and from the 
lemma lA~2l we have: 



Proof Of Theorem [TJ 

Substitute the control laws ([7]) into the equation: 

m 

rh(t) = O a m{t) + E UiO^mit) + Dm(t) + g. (A5) 

i=l 

Let m if), m 2 (t) are two solutions of the equation (|A5|) 
with the initial values to be mj, respectively, then 
rn 1 (t) — m (t) satisfies the following time- variant linear 
equation: 

rn 

±{t) = (D + O + E UiO?)x(t) (A6) 



with the initial value x(to) = — TOq. 

Let Q = I, rj = fi = 1, u = 2d mm , where -d mm < 
(from the assumption (H2)) is the maximal eigenvalue of 
D. From lemma [A~T1 it can be verified that: 

|m x (t) - m 2 (t)\ < e - 2d "' i "< t - to Vo - m 2 \, 

which means wn}(t) — m 2 (t) — > when t tends to infinity. 

From the above analysis, to show that the solution of 
equation © tends to §8§ , it is sufficient to prove that ([5]) 
satisfies the equation (|A5[) . which is equivalent to 



E "iO? e^-^m^ + De° o{t ^m^ + g = 0, 



(O?m 0O (t),-..,Qg,m 0o (t)) 



e -0 (*-to)(9P e Oo(t-to) m o0 5 . . . j 

-O (t-i )Qp e° o( *~ to) mg°) e -0"(*-*o)^ 

(771 
E(e"° !1(( " to) )uO^<,- 
fe=l 



fc=l 



= e Oo(*-*o)(Ojm~,... ,OP l mg )(e- o 1 ( t - t «>)) r e- "( t - t 



e Oo( f - to )^ 60 P TOo 



Therefore, the equation (|A7|) is reduced to 



e Oo( t ~to) g e . P TO oc + Dm oc + g j = Q) 



which is equivalent to 
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